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Abstract. 

Optical interferometry is a powerful tool to investigate the close environment of 
AGB stars. With a spatial resolution of a few milli-arcseconds, it is even possible to 
image directly the surface of angularly large objects. This is of special interest for Mira 
stars and red supergiants for which the dust-wind is initiated from or very close to the 
photosphere by an interplay between pulsation and convection. Based on two-epoch 
interferometric observations of the Mira star X Hya, we present how the variation of 
the angular size with wavelength challenges pulsation models and how reconstructed 
images can reveal the evolution of the object shape and of its asymmetric structures. 


1. Introduction 


The mass loss that determines the fate of AGB stars is triggered by convection and 
pulsation phenomena. The subsequent formation and transport of material from the 
surface to the interstellar medium still remain to be discovered. The very high resolu¬ 
tion that interferometric imaging delivers allows to observe the surface, the molecular 
layers as well as the dusty atmosphere. In particular, asymmetries in the brightness 
distribution, e.g. spots on the surface or clumps in the atmosphere, are the signatures of 
processes such as convection or pulsation. The key to pin down their nature and origin 
is to derive their spectral dependence and their evolution timescales. These properties 
are guaranteed to bring decisive insights on the mass loss processes at play in these 
stars. With this aim, we here report on interferometric observations of the mira star X 
Hya observed with AMBER. 
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2. Spectro-interferometric observations 


X Hya is a M7e Mira star w ith a 301 day pulsation period loc ated at a distance o f ~ 440 
pc ( Whitelock et al.1 12008 ). We observed it with AMBER dPetrov et al.1 2007) at two 
epochs corresponding to two phases of its pulsation cycle: ® = 0.0 and ® = 0.2. A 
multi-wavelength observation is required to distinguish the counterpart of the molecu¬ 
lar layer compared to the stellar photosphere. We thus performed J-H-K observations 
(covering from 1.17 gm to 2.37 gm) at low spectral resolution (R ~ 30). 

The datasets we obtained at each epoch are rich and complex. Even at short 
baselines, squared visibilities ( V 2 ) show a significant departure from a simple centro- 
symmetric geometry such as an uniform disk model (Fig. [l}. In order to understand the 
spectral and morphological signatures on the V 2 data, we first divided the data in six 
spectral bands as shown on Tab. [T] This division is based on an analysis of the visibility 
variation with the wavelength and is coherent with a previous spectro-interferometric 
analysis of X Hya dWittkowski et alfeoi ll) . 

We first performed a uniform disk (UD) model to get a first estimate of the an¬ 
gular size of the object in the different bands and at the different position angles (PAs) 
we observed. We only considered squared visibility data whose spatial frequencies be¬ 
longed to the first lobe in order to avoid any contamination from limb darkening or 
structures smaller than the stellar diameter. An average of the measurements over the 
PAs is shown in Table [0 
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Figure 1. Squared visibilities data at Epoch A (left) and Epoch B (right). Colors 
indicate the 6 different spectral bands. The solid lines represent a 6 mas (for Epoch 
A) and a 7 mas (for Epoch B) uniform disks at 1.75 microns (average wavelength of 
the covered spectral interval). 


2.1. Radiative transfer modeling 

To better understand the spectral dependence of this complex dataset, we used a grid of 
dynamic model atmosphere series computed with the P/M and CODEX codes. These 
models develop a self-excited pulsation mechanism and sample opacities over 4300 
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Wavelength range 

Contribution 

< Pud 
<t> - 0.0 

4>ud 

O = 0.2 

1.1 - 1.4/on 

Continuum (1) 

5.94 ± 0.15 mas 

7.41 ± 0.40 mas 

1.4 - 1.5 /mi 

H 2 0(1) 

6.90 ±0.13 mas 

7.85 ± 0.49 mas 

1.5 - 1.7 /mi 

CO(l) 

6.03 ± 0.59 mas 

7.19 ± 0.85 mas 

1.7 - 2.1 /mi 

H 2 0 (2) 

6.52 ± 0.23 mas 

8.24 ± 0.86 mas 

2.1 - 2.25 /mi 

Continuum (2) 

6.20 ± 0.35 mas 

7.02 + 0.86 mas 

2.25 - 2.4/mi 

CO (2) 

6.80 ± 0.42 mas 

7.54 ± 0.79 mas 


Table 1. Division of the covered spectral range in 6 spectral bands and uniform 
disk diameters for the two observing epochs. 



Spatial frequency (1/arcsec) 


Figure 2. Modeling of the squared visibility data (in black) with CODEX/R52 
(red), P (in blue) and UD (cyan) models. The Rosseland diameter indicated in the 
legend was derived from the best-fit of these models. 


wavelengths running from 200 nm to 50 //m, with 1 nm sampling between 200 nm and 
3 /mi. For further details scc llreland et al.l ( 2008 1 1201T ). 

For a set of parameters, these models produce radial profiles of the intensity across 
the chosen spectral band that are turned into V 2 by a Fourier transform. These latter 
were finally fitted to our data by adjusting the Rosseland diameter and the pulsation 
phase (for CODEX models only). Figure |2] shows V 2 -fitted curves for one baseline. 
At the difference of the simple geometric models, these physical models that take into 
account Hot) and CO opacities, do reproduce the variation of the V 2 with the spatial 
frequency fairly well. The results of the model parameters fitting arc summarised in 
Tab. [2] 


3. Image reconstruction 

In the previous sections, we only took into account the V 2 data which means the object 
was thought to be centro-symmetric. Flowever, because AMBER recombines the beams 
from three telescopes, we do have a phase information called closure phase that is a 
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Models 

Epoch A, <J> = 0.0 

Epoch B, ® = 0.2 

Best P model 

P27n 

P14n 

Gross (mas) 

4.95 

7.38 

xLi 

11.5 

62.1 

Best M model 

M16n 

M12 

Gross (mas) 

5.00 

6.58 

xLi 

14.3 

63.2 

Best R52 model 

388040 

386280 

Gross (mas) 

4.02 

4.54 

y 2 

"■ red 

8.1 

44.3 

Best o54 model 

O54/262160 

O54/287940 

6 Ross (mas ) 

4.76 

5.86 

X 2 red 

7.4 

41.2 

®CODEX 

0.7 

0.5 


Table 2. Results of the fitting of radiative transfer models. 


measurement of the asymmetry of the object. The closure phase signals exhibited a 
significant departure from the null value in all spectral bands. 

In order to reconstruct the spatial intensity distribution out of such a complex 
dataset of V 2 and closure phases, it is mandatory to use image reconstruction algo¬ 
rithms such as MIRA dThiebauti 120081) . However, given our limited UV coverage in 
individual bands, we can’t aim at reconstructing a model-independent image. Nev¬ 
ertheless, we carried out an image reconstruction approach based on our parametric 
modeling results. We therefore used ID intensity profiles from the radiative transfer 
models that we turned into a prior image. The reconstructed image was obtained by 
adjusting the image to the data and by regularizing towards the prior image. 

Several regularization techniques were used. The best results were obtained with a 
quadratic regularization towards an a priori object image derived from the R52/388040 
model intensity profile. All the reduced^ 2 values associated with the best reconstructed 
images are lower than 1. Figure [3] shows images for two spectral bands at the two 
phases. 

In agreement with our previous results, the overall size of the object increased 
between the two epochs in all spectral bands. X Hya also presents a more complex 
morphology at the second phase and especially in the continuum and in the H 2 O bands. 
In order to push further the interpretation of these observations, we plan to compare 
these images with full 3D hydro-radiative transfer models that also include self-excited 
pulsating atmospheres. 


4. Conclusion 

We report on low-resolution spectro-interferometric AMBER observations of the mira 
star X Hya. After a first simple geometric analysis, we performed a model-fitting with 
the P/M and CODEX atmosphere series that were able to reproduce the chromatic de¬ 
pendence of the squared visibilities. In order to adjust the non-zero closure phase sig¬ 
nals that we observed in all spectral bands, we used the reconstruction algorithm MIRA 
to obtain images at the two phases of the pulsation cycle. We observed a higher degree 
of asymmetry at the second epoch (® = 0.2) in the continuum and in the H 2 O bands. 
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Figure 3. Reconstructed images with MIRA for spectral bands ’’Continuum 1” 
and ’’H 2 O (2)” at Epoch A (left column) and Epoch B (right column). The inter¬ 
ferometric beam size is displayed in the bottom-right corner of each image. White 
contours mark areas of 25%, 75% and 95% of the maximum flux of the image. 


We will compare these reconstructed images with 3D hydro-radiative transfer models 
of pulsating atmospheres. 
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